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The heavy gravitino in the minimal supergravity (mSUGRA) models is likely to be the lightest 
supersymmetric particle (LSP). Produced from the late decays of the metastable Weakly Interacting 
Massive Particles (WIMPs) such as the lightest neutralinos, the stable gravitinos can be plausible 
candidates for the cold dark matter in the universe. Such gravitino dark matter can naturally 
evade the current detection experiments due to its superweak couplings. However, this scenario 
must be subjected to the constraints from the Big Bang nucleosynthesis (BBN) predictions for 
light element abundances as well as the Wilkinson Microwave Anisotropy Probe (WMAP) data 
for the relic density. Assuming the popular case in which the lightest neutralino is the next-to- 
lightest supersymmetric particle (NLSP), we find that requiring BBN predictions for light element 
Oh' abundances to agree with the WMAP data can impose upper and lower mass bounds on both the 

gravitino LSP and the neutralino NLSP. A scan over the mSUGRA parameter space, subjected 
to the BBN constraints, the WMAP data and the b — > S7 bounds, shows that the low tan/3 
' 40) region as well as the region accessible at CERN Large Hadron Collider (LHC) will be severely 

£N) , constrained. Such stringent constraints on the parameter space might be instructive for testing this 

scenario in future collider experiments. 
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•£} I. INTRODUCTION 
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The nature of the dark matter is one of the mysteries in today's physical science. It has been intensively explored 
both theoretically and experimentally. Studies showed that the cosmic dark matter is plausibly composed of non- 
baryonic Weakly Interacting Massive Particles (WIMPs) [1]. While the Standard Model of particle physics cannot 
provide a candidate for the dark matter WIMP, the popular supersymmetric theory with R-parity conservation can 
Oh. provide a good candidate, i.e., the lightest supersymmetric particle (LSP). So far the widely studied scenario is that 
the lightest neutralino is assumed to be the LSP. However, despite of the overwhelming popularity of this scenario, 
other possibilities should not be ignored due to the following reasons. On the one hand, the success of such neutralino 
dark matter scenario may be spoiled by the problems caused by gravitino in the reheating era [2]. On the other hand, 
the neutralino dark matter scenario has not yet been confirmed by current experiments [3]. 

One possible scenario other than neutralino LSP is that the gravitino is assumed to be the LSP. Such gravitino 
LSPs can form warm or cold cosmic dark matter, depending on the gravitino mass: 

(i) In some low-energy SUSY breaking models, like the gauge mediated SUSY breaking (GMSB) models, the 
gravitino can be as light as KeV, much lighter than other supersymmetric particles. It can thus form the warm 
dark matter. Note that the recent WMAP data imposed severe constraints on the dark matter type. As analyzed 
in [4], while a very tiny component of dark matter can be the hot neutrinos, the warm dark matter is ruled out 
due to the detected early re-ionization of the universe at a redshift z w 0.20. Therefore, the scenario of warm 
dark matter gravitino is not favored by recent observation. 

(ii) In the popular mSUGRA models, the gravitino mass is unspecified and only known to be of the weak-scale. 
Such heavy gravitino is possibly the LSP and can form the cold dark matter in the universe. In contrast to 
the highly constrained scenarios, in which the gravitino is produced as a thermal relic [5] or produced during 
reheating [6], a new scenario, assuming the gravitino to be produced from the late decays of the thermal relic 
WIMPs, was recently proposed in [7]. 

In this article we focus our attention on this new gravitino dark matter scenario. Since the heavy gravitinos couple 
gravitionally, they are naturally the so-called Superweakly Interacting Massive Particles (SuperWIMPs). As a plausi- 
ble candidate for the cold dark matter in the universe [7], the gravitino Super WIMP can naturally evade the current 
dark matter detection experiments due to its superweak couplings. However, this scenario must be subjected to the 
constraints from the Big Bang nucleosynthesis (BBN) as well as the WMAP data [8] : 
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• The late decays of WIMPs (like neutralinos) into gravitino SuperWIMPs will release electromagnetic (EM) 
and hadronic energy. Such energy release will alter the BBN predictions for light element abundances [9,10]. 
Requiring the resulted predictions for light element abundances to agree with the measured values will impose 
strong constraints on the gravitino dark matter scenario. 

• WMAP precisely measured many quantities, especially the total matter density and the baryon density, 

n m h 2 = 0.1351°°°* , n b h 2 = 0.0224t°;°°g . (1.1) 

From such results we can deduce the 2a range for the cold dark matter density 

VcDMh 2 = 0.1126±g:8i|l , (1.2) 



which is dramatically more accurate than previous results and agrees quite well with other approaches. Such 
precise measurements will impose strong constraints on the parameter space of gravitino dark matter scenario. 

The aim of this article is to examine these constraints on the parameter space of mSUGRA in this new gravitino dark 
matter scenario. Assuming the popular case that the lightest neutralino is the NLSP, we will examine the constraints 
on the mSUGRA parameter space from the BBN light element abundances, the WMAP data of relic density as well 
as the b — > branching ratio data. 



II. CONSTRAINTS ON THE PARAMETER SPACE 

In the new scenario [7], the NLSPs freeze out with a thermal relic density Qnlsp and then decay to the gravitino 
at time 10 4 <~ 10 8 s. Thus the relic density of gravitino dark matter is obtained by 

The late decays of the NLSPs will release energy which will alter the light clement abundances. In fact, the later 
injection of high energy photon with the stopping energy inversely proportional to the temperature by scattering off 
the background photon will dissociate the existing light elements. If carefully chosen, such injecting EM energy can 
destruct the light element abundances to proper values. It is well known that the BBN predictions for light element 
abundances are quite successful for most light elements. Yet the BBN predictions for 4 Hc and 7 Li (especially 7 Li) 
seemingly do not agree with the WMAP data [11]. Requiring such discrepancy be settled by the energy released 
from the late-decaying NLSPs and, at the same time, requiring such energy release not to spoil the successful BBN 
predictions for other light elements, the constraints on the NLSP lifetime (in seconds) and EM energy release can be 
obtained [9]: 

1.5 x 10 6 s < t < 4 x 10 6 s, (2.2) 
0.8 x 10~ 9 GeV < Cem < 1.5 x 10~ 9 GeV, (2.3) 

where 

Cem = eem-Bem^nlsp (^-4) 

is the emitted EM energy density, with 6em being the initial EM energy release from the decay, -Bem being the 
branching fraction of the decay into EM components, and Ynlsp = «nlsp/% being the NLSP number density 
normalized to the BG photon number density. In the derivation of the above bounds, the ratio r\ = ns/n^ is fixed to 
be 6 x 10~ 10 , and the bounds on the primordial abundances of light elements are taken to be [9] 

1.3 x 10~ 5 < D/H < 5.3 x 10" 5 , (2.5) 

0.227 < Y p < 0.249 , (2.6) 

9.0 x 10" 11 < 7 Li/H < 2.8 x 10~ 10 , (2.7) 

6 Li/ 7 Li <, 0.07 , (2.8) 

6 Li/H <, 2 x 10" 11 , (2.9) 

where Y p denotes 4 He abundance. 
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The dominant decay of neutralino (%) NLSP into gravitino (G) LSP is through \ — * 7^ with a rate given by 

\N X1 cos9 w + N 12 sin9 w \ 2 m| 



r(x - iG) = 
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where is the gravitino LSP mass, is the neutralino NLSP mass, M* = l/-\/87rGAr ~ 2.4 x 10 18 GeV is the 
reduced Plank scale, Qw is the weak mixing angle, and Nij denotes the matrix element projecting the i-th neutralino 
into Bino (j = 1), Wino (j = 2) and Higgsinos (j = 3,4). For such a decay we have 
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We would like to make some clarifications about our numerical calculations: 



(1) In the calculation of the thermal relic density of the NLSPs, we considered the general mixing of neutralinos 

and used the package Microomega [13], which includes all tree level contributions to the scattering amplitudes 
l 

(2) In the scan over the mSUGRA parameter space, we used SuSpect2.0 [16] to obtain the sparticle masses in 
mSUGRA models, which includes one-loop corrections to sparticle masses and two-loop corrections to Higgs 
masses. 

(3) When deriving the constraints on the mSUGRA parameter space, we consider the BBN constraints in 
Eqs.(2. 2,2.3), the WMAP data of relic dark matter density in Eq.(1.2) as well as the b — > sj bound [17]: 



2.16 x 1CT 4 < BF(b -» sj) < 4.34 x 10" 



(2.13) 



Furthermore, we fixed the common trilincar coupling Aq = since our results are not sensitive to it. 

(4) In our calculation we also considered the muon anomalous magnetic moment a^. Since so far much theoretical 
uncertainty exists in predictions 2 , we did not use a M to set constraints and, instead, we only gave the range 
of the mSUGRA contributions corresponding to each case. This will be useful when a M uncertainty is further 
reduced in the future. 

The bounds on the masses of the neutralino NLSP and gravitino LSP are shown in Fig.l. One sees that BBN results 
give rather stringent upper and lower bounds 

430 GeV < M LS p < 600 GeV, 630 GeV < M NLS p < 840 GeV. (2.14) 

We found that the neutralino NLSP is quite bino-like (Nu ^ 0.99) in the allowed regions in Fig.l 

The allowed regions in the plane of m versus m 1 / 2 are shown in Fig. 2 for plus sign of \i and Fig. 3 for minus sign of 
[i. From Fig. 2 we see that for plus sign of fi much of the parameter space with tan/? ^ 40 is ruled out. For example, 
for tan/3 = 10, only a narrow strip survives, i.e., 225 GeV <, m <y 300 GeV and 1020 GeV <> m 1/2 <> 1270 GeV (the 
corresponding contribution to is 1.38 ^ Aa M x 10 10 2.05). When tan/3 increase to 50, the allowed region gets 
quite large, i.e., 740 GeV ^ mo ^ 1730 GeV and 1060 GeV ^ ^ 1720 GeV (the corresponding contribution to 
a M is 2.18 <, Aa M x 10 10 £ 6.32). 

For minus sign of /i, similar results are obtained. We see from Fig. 3 that for tan/3 30, the constraints are quite 
stringent. For example, for tan/3 = 10 the constraints are 220 GeV m ^ 305 GeV and 1020 GeV ^ m 1 / 2 SS 
1280 GeV (the corresponding contribution to is —6.07 ^ Aa p x 10 10 ^ —3.7). When tan/3 increases to 40, 
the constraints are weakened to 1060 GeV ^ m ^ 1530 GeV, 1290 GeV ^ mi/ 2 & 1795 GeV (the corresponding 
contribution to is —5.45 ^ Aa^ x 10 10 ^ -1.90). 



1 In the present version of Microomega, the package FeynHiggs [14] is used to calculate the Higgs masses and Hdecay [15] is 
used to include important QCD corrections to the Higgs decays 

2 For example, one group [18] gives 11.5 < Aa M x 10 10 < 60.7 or -16.7 < Aa„ x 10 10 < 49.1, depending on the calculation 
approaches. Another group [19] gives <5a M x 10 10 = 12.41^3 by using r decay data to determine vacuum polarization. 
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FIG. 1. The region (shaded) allowed by BBN to account for 7 Li and 4 He abundance in the scenario of gravitino LSP and 
neutralino NLSP. 
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FIG. 2. The regions (shaded) allowed by the WMAP data, BBN constraints as well as 6 — > s*y data. From left to right 
corresponds to tan /3 = 10, 20, 30, 40, 50, respectively. The sign of fi is assumed to be plus. 
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FIG. 3. Same as Fig. 2, but for minus sign of fx. 

It is interesting to note that if Aa^ is required to be positive (negative), then the allowed regions in /.i < (fi > 0) 
case will be completely ruled out. Therefore, a further scrutiny of a M will provide meaningful constraints on the 
parameter space of mSUGRA. 

It should be pointed out that the above stringent constraints in Figs. 2 and 3 were obtained under the requirement 
that the gravitinos from the late decays of the NLSP make up all the dark matter. If other particles like axions are 
also a component of dark matter, then the constraints will be changed. Moreover, if the relic gravitino dark matter 
are partly produced during reheating 3 , the constraints will also be changed. In order to forbid gravitino (with mass 
range in Eq.(16)) production during reheating, the reheating temperature Tr must be lower than a certain value. 
Following the analyses in [6], we evaluated such a limit and found Tr 10 9 GeV. 

Note that in our analysis we only considered the consequence of EM energy release from the late decays of WIMPs 
and required such EM energy release to settle the discrepancy between the BBN predictions and the observed values 
for light element abundances. However, as can be seen in the allowed Mnlsp/Mlsp region of Fig. 1, the decay 
NLSP — > G + Z is also possible and thus will cause hadronic energy release from the subsequent Z decays. Such 
hadronic energy release can be quite dangerous since it can alter the BBN predictions for the light element abundances. 

In Ref. [21] the authors studied the bounds from the hadronic energy release and found that such hadronic bounds 
may be stringent 4 due to the overproduction of D. However, as pointed out in Ref. [21], for the best fit region 
with lifetime r between 10 6 s and 10 7 s, the D overproduction from hadronic energy is possibly canceled by the 
overdestruction from EM energy and thus the hadronic bounds become less stringent and subject to large uncertainty. 
In such a region the EM energy effects are dominant [21]. Since in our study we focused on the best fit region, we did 
not consider the hadronic effects. 



III. CONCLUSIONS 

We examined the constraints on the newly proposed dark matter scenario, in which the gravitino is assumed to be 
the LSP and produced from the late decays of metastable NLSP. Although such gravitino dark matter can naturally 
evade the current detection experiments due to its superweak couplings, we found, however, that this scenario is 
subjected to stringent constraints from the BBN predictions for light element abundances as well as the WMAP data 



3 The gravitinos can also be produced as a thermal relic at very early times. However, in the context of inflation, the universe 
inflated between that early time and now, which would dilute any gravitino thermal relic density. 

4 In their study they used the bounds 2.4 x 10 -5 < D/H < 3.2 x 10~ 5 , which are more aggressive than the bounds of Eq.(2.5) 
used by us to obtain the best fit region. 
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for the relic density. Assuming the popular case that the lightest neutralino is the NLSP in mSUGRA models, we 
found that the low tan (3 40) region as well as the region accessible at the LHC are severely constrained. 

The popular mSUGRA models will be explored in future colliders like the LHC. For this purpose, it is important to 
know which part of the parameter space is viable and thus should be primarily explored. In this regard, the stringent 
cosmological constraints on the mSUGRA parameter space obtained in this work will be useful. Especially, if the 
ongoing and planned dark matter detection experiments fail to find any dark matter signal, it will imply that the 
dark matter interactions are too weak and thus the gravitino dark matter scenario will be favored. Then, to test this 
scenario at colliders, the stringent cosmological constraints on the parameter space will be quite instructive. This 
would serve as a good example that the studies in cosmology and astrophysics can shed some light on particle collider 
physics. 

On the other hand, the LHC could explore mSUGRA parameter space up to m 1 / 2 ~ 1400 GeV (700 GeV) for 
small (large) values of mo, assuming 100 of integrated luminosity [20]. If the LHC results finally restrain the 

parameter space to one of the regions obtained in this work, then it implies that the gravitino dark matter scenario is 
favored. In this sence, the studies in particle physics can provide some insights in the understanding of dark matter 
in cosmology. 

We address that our study in this work is just illustrative instead of exhaustive. We assumed the popular case that 
the lightest neutralino is the NLSP in mSUGRA models. Actually, other super particles, like tau-slepton, are also 
likely to be the NLSP in mSUGRA models. If tau-slepton is assumed to be the NLSP, there are some theoretical 
uncertainties in its decay modes and the corresponding energy release. 

Note added: While we are preparing this manuscript, some other preprints [21] appeared, where the constraints on 
the gravitino dark matter scenario are studied. We found that the studies in [21] are quite exhaustively, where the 
scenarios of neutralino NLSP, stau NLSP as well as sneutrino NLSP are all considered. Compared with the studies 
in [21], the characteristic of our study is that we performed a scan over the mSUGRA parameter space and presented 
the allowed regions in terms of original mSUGRA parameters. In addition, the BBN constraints on the EM energy 
release are more stringent in our study since we required such energy release settle the discrepancy between the BBN 
prediction and the WMAP data for 7 Li abundance. 
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